Parameter study of the diamagnetic relativistic pulse accelerator (DRPA) in slab 
geometry I: Dependence on initial frequency ratio and slab width 
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Two-and-a-half-dimensional particle-in-cell plasma simulations are used to study the particle en- 
ergization in expanding magnetized electron-positron plasmas with slab geometry. When the mag- 
netized relativistic plasma with high temperature (initial electron and positron temperature are 
ksTe — ksTp = 5MeV) is expanding into a vacuum, the electromagnetic (EM) pulse with large 
amplitude is formed and the surface plasma particles are efficiently accelerated in the forward di- 
rection owing to the energy conversion from the EM field to the plasma particles. We find that the 
behavior of the DRPA (Diamagnetic Relativistic Pulse Accelerator) depends strongly on the ratio 
of the electron plasma frequency to the cyclotron frequency Upe/fle and the initial plasma thickness. 
In the high Upe/Qe case, the EM pulse is rapidly damped and the plasma diffuses uniformly without 
forming density peaks because the initial thermal energy of the plasma is much larger than the 
field energy. On the contrary, in the low ujpe/^le case, the field energy becomes large enough to 
energize all the plasma particles, which are confined in the EM pulse and efficiently accelerated to 
ultrarelativistic energies. We also find that a thicker initial plasma increases the maximum energy 
of the accelerated particles. 

PACS numbers; 
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I. INTRODUCTION 

In astrophysical, solar and space plasmas, the accel- 
eration of high energy particles is an intriguing and un- 
solved problem. Magnetic reconnection, wave turbulence 
and coUisionless shocks are often invoked for the acceler- 
ation of plasma particles, from cosmic rays, solar flares 
to gamma-ray bursts,!] and astrophysical jets. 

A new mechanism of particle energization by freely ex- 
panding, hot, strongly magnetized coUisionless plasmas 
was recently discovered by Liang et al. using two-and-a- 
half dimensional (2-1/2D) particle-in-cell (PIC) simula- 
tions 2]. This mechanism, called the diamagnetic rela- 
tivistic pulse accelerator (DRPA) , is able to convert most 
of the initial magnetic energy into the ultra-relativistic 
directed kinetic energy of a small fraction of the sur- 
face particles. Moreover, Liang and Nishimura|3| have 
found that the late-time plasma pulse of the DRPA repro- 
duces many unique features of cosmic gamma-ray bursts. 
Hence the DRPA may shed important new light on the 
unsolved problem of gamma-ray bursts. When a hot, 
strongly magnetized (/3=thermal pressure/plasma pres- 
sure < 1) coUisionless plasma expands into a vacuum 
or low-density surrounding, an electromagnetic pulse is 
formed at the expansion surface, which traps and accel- 
erates a small fraction of surface particles via the ponder- 
motive forced, H, and the magnetic energy is efficiently 
converted into directed particle energy . 

In our previous Letters 1113, the DRPA was studied 
by performing 2-1/2D PIC simulations for both electron- 
positron (mass ratio is nii/me — 1) and electron- ion 
{mi/rrie — 100) plasmas with a specific set of initial 
parameters, and only carried out to lOOOil^T^ (where 
rie = eBo/rrie is the electron cyclotron frequency de- 



fined by the initial magnetic field). In this paper, we 
will study the parameter dependence and the long-term 
evolution of the DRPA of electron-positron plasmas by 
means of a series of large-scale PIC simulations. 

The outline of this paper is as follows. Section ^ de- 
scribes the simulation model and the simulation param- 
eters. The benchmark results of Liang et hLQ] are re- 
viewed and summarized in Section IIII Al for comparison 
with other cases. The long-term evolution of the DRPA is 
described in Sec. IIII Bl The initial parameter dependence 
of DRPA is discussed in Sees. ImUl and UlIDl Section ITTI 
gives a summary of our results. 



II. SIMULATION MODEL 

In our particle simulation codeQ, we use a 2-1/2D ex- 
plicit simulation scheme based on the PIC method for 
time advancing of plasma particles. In this method, spa- 
tial grids are introduced to calculate the field quantities, 
and the grid separations are uniform. Ax — Az — Ae, 
where Ae is the electron Debye length defined by c/tOpe 
(c is the speed of light and the electron plasma frequency 
is ujpe = -y/ e^rio/eoTTie; e is elementary charge, uq is the 
initial electron density, and eo is the dielectric constant 
of vacuum). In our simulations, the thermal speed of 
electrons is almost equal to c because the plasma tem- 
perature is ksTe = ksTp = 5MeV, where the subscripts 
e and p refer to electrons and positrons. The simulation 
domain on the x ~ z plane is —Lx/2 < x < L^/2 and 
~L,/2 <z< 

We treat plasma particles as superparticles, with 20 su- 
perparticles per cell (if the superparticles are spread uni- 
formly on the entire simulation domain) for each compo- 
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nent (electrons and positrons) . Most of the particle den- 
sities, however, concentrate on the expanding surface in 
our study and we are mainly interested in particle behav- 
iors near the surface. Thus, this number of superparticles 
is enough to accurately resolve the physics of surface par- 
ticle acceleration. To prevent large violations of Gauss's 
law due to numerical noises, the Marder's method for the 
electric field correction 7J is adopted at every time step 
in the code. 

We use a doubly periodic system in x and z direc- 
tions, and the system length is Lx = 3840 Aa; (Run A), 
3840Aa; (Run B), 1920Ax (Run C), and 7680Aa; (run D) 
and Lz = 10 Az (for all four Runs). Initially, the elec- 
tron and positron distributions are assumed to be rela- 
tivistic Maxwellian with spatially uniform temperature, 
fc^Tg = ksTp — bMeV . The spatial distribution of ini- 
tial plasmas is assumed to be a uniform slab with the 
lengths I2A2; X 10Az,(runs A,B,C) and 60Ax x lOAz 
(run D) and the plasmas are located in the center of the 
system. The background magnetic field = (0, i3o,0) 
initially exists only inside the plasma. 

We will first compare three simulation results by 
changing the frequency ratio ujpf,/^e- The values are 
ujpe/^e = 0.104 (Run A), 1.0 (Run B), and 0.01 (Run 
C). At the end we will also discuss the result of changing 
the initial slab thickness (Run D). 



III. RESULTS 
A. Diamagnetic relativistic pulse accelerator 

First, we illustrate the mechanism of the diamagnetic 
relativistic pulse accelerator (DRPA) by reviewing the 
early {tVle < 10'^) results of the benchmark Run A (Liang 
et ai, 20030)- In Run A, the magnetic field energy is ini- 
tially almost equal to the thermal energy of plasma par- 
ticles, i.e., electron and positron beta is f3e = Pp = 0.551 
(/3 = 2^QnokBT/ Bq, where /^o is the magnetic perme- 
ability of vacuum). 

In Figure nia), the system-integrated energies of the 
magnetic field (line 1), electric field (line 2), electrons 
(line 3), and positrons (line 4) are shown as functions of 
time. The initial static magnetic field rapidly transforms 
into two oppositely propagating electromagnetic pulses, 
which propagate towards the vacuum region. The plasma 
simultaneously begins to expand. We see the rapid de- 
crease of the magnetic energy and the increase of the 
electric field energy due to the generation of the electro- 
magnetic (EM) pulse during the period shown as Phase 1 
in Fig. ^a). In Phase 2, both the magnetic and electric 
field energies begin to decrease and they are converted 
into the electron and positron kinetic energies. Also, the 
EM pulse with the amplitude cBy ~ is formed at the 
expansion surface during this period. In Phase 3, the 
energy conversion from the field energies to the particle 
kinetic energies still continues although the conversion 
rate becomes slower compared with that in Phase 2. 
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FIG. 1: System-integrated energy in the magnetic field (line 
1), electric field (line 2), electrons (line 3), and positrons (line 
4) as functions of time for (a) Run A, (b) B, (c) C, and (d) D. 
Each energy component is normalized by the initial magnetic 
field energy in each figure. 



Figure El illustrates the spatial profiles of the (a) mag- 
netic field By, (b) electric field E^, (c) current density 
Jz, and (d) electron density rig at tQe — 0, 80, 500, 
and 1000 for Run A. The magnetic field and the elec- 
tron density for a; < are identical, on the other hand, 
the electric field and the current density become antisym- 
metric about X = 0. The values of Ez and Jz are zero 
at tQe = 0. Only the profiles near the expansion surface 
are drawn at tile = 500 and 1000 to avoid making the 
figure cluttery. At tfle = 80, the EM pulse is formed 
and begin to expand with plasmas into the vacuum. The 
gradient drift in the z-direction induces a current den- 
sity. At trie = 500, the EM pulse has been formed at the 
expansion surface. Most of the plasma particles follow 
the EM pulse propagation and the plasma density has 
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FIG. 2: Spatial profiles of the (a) magnetic field By, (b) elec- 
tric field Ez, (c) current density Jz, and (d) electron density 
rie as functions of a; at 2 = for Run A. The quantities only 
in the a; > region are shown at tQ,^ — 0, 80, 500, and 1000 
(reproduced from Ref.[2]). 
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FIG. 3: 2-D Electron distribution on the simulation plane 
for Run A at (a) tQ.e = 0, (b) 500, and (c) 1000. 



a peak corresponding to particles trapped by the pon- 
dermotive force of the pulse. The amplitude of the EM 
pulse is gradually damped as EM energy is transformed 
into particle energy. 

The X — z configuration of electrons is plotted at three 
different times in Figure |31 The electron-positron slab 
with the internal magnetic field expands rapidly into the 
vacuum region. The behavior of positrons is identical 
with that of electrons. The expansion speed of the plasma 
surface in the x direction is almost the speed of light, 
Vx — c. Several simulation studies of non-relativistic 
plasma expansion across an initially uniform magnetic 
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FIG. 4: Comparison of phase plots of electrons in the (a) 
X — Px and (b) x — Pz spaces for Run A at tO.^ — 0, 80, 500, 
and 1000. 



field have been carried out[S|,|2l,lli]|,llH, and the lower hy- 
brid drift instability has been found to deform the struc- 
ture of the expanding plasma surface. As shown in Fig.O 
there appears no instability at the plasma surface in our 
result. There are two main differences between our sim- 
ulation and the previous studies |S 1^ ITol ITU : (a) in their 
studies, an initial magnetic field exists uniformly all over 
the simulation plane and the plasma expands (or diffuses) 
across the field, (b) the electron-ion plasma was consid- 
ered in their studies. Thus, the charge separation due to 
the large mass ratio causes the different drift speed be- 
tween electrons and ions, which makes the lower hybrid 
drift wave unstable. These differences may explain the 
absence of instability in our simulation. In addition, the 
effective Larmor radius of relativistic plasmas given by 
7(c/f2e) is larger than our system length in the z direc- 
tion at late times. Thus, any effect of the Larmor radius 
play no important role in our simulations, and we cannot 
see any effect caused by the Larmor motion of relativistic 
plasmas. 

Figure 0] shows phase plots of electrons at four dif- 
ferent times in the space (a) x — and (b) x — Pz 
for Run A. It is found that most of the electron pop- 
ulation moves to the expansion surface and the momen- 
tum Px of the surface electrons increases as time elapses. 
The averaged momentum reaches < Px >— lOOrTieC and 
the maximum momentum of energetic particles reaches 
Px,max ~ 200r7ieC at iSlg = 1000. On the other hand, 
the momentum P^ indicates no dramatic increase even 
at late times. This means that the plasma particles are 
effectively accelerated into the x direction rather than 
the z direction. 
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FIG. 5: Spatial profiles of the (a) magnetic field By, (b) elec- 
tric field Ey, (c) current density Jz, and (d) electron density 
rie as functions of a; at z = 0, and (e) phase plots of electrons 
in the x — space for Run A. The quantities only in the 
a; > region are shown at tQ.e = 6000, 10000, 14000, and 
18000. 



B. Long-term evolution of the DRPA 

Let us focus on the long-term evolution of the DRPA 
in Run A. Figure O shows the spatial profiles of the (a) 
magnetic field By, (b) electric field Ez, (c) current density 
Jz, (d) electron density Ue, and (e) the phase plots x — P^ 
of electrons for Run A at = 6000, 10000, 14000, and 
18000. Only the portion near the expansion front at each 
time has been displayed in each panel. In Figs.|3a) and 
Ob), the wave packet is dispersively spreading and the 
EM pulse amplitude gradually decreases as time elapses. 
The reduction of the pulse amplitude corresponds to the 
conversion of the EM energy into the kinetic energy of the 
expanding plasma. In Fig.jSfe), we recognize the electron 
momentum Px near the first pulse of the wave packet 
dramatically increases as time goes on. The maximum 
value of Px reaches more than SOOrTigC at tfie = 18000. 

Another significant feature of the late time behavior 
in Run A is the bifurcation of plasma density. In the 
initial phase of the DRPA, the electron density near the 
wave packet has only one peak as shown in Fig. [2Jd) . At 
late times in Fig.l^d), the peak has been split into sev- 
eral peaks. The bifurcation of the plasma density profile 
results from the development of new transverse current 
and reverse currents behind the front, which forms new 
traps for particles in the pondermotive potential. 

FigureEfa) shows the momentum distributions of elec- 
trons which are located in the vicinity of the EM pulse 
at = 1000, 6000, and 10000. The momentum peak at 
which the distribution takes a maximum value tends to 
increase as time elapses, i.e., Px — 4:5meC at t^le = 1000, 
Px ~ 60TOeC at ifie = 6000, and Px ~ 80TOeC at tfle = 
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FIG. 6: Electron momentum distributions for (a) Run A, (b) 
B, (c) C, and (d) D. Each distribution consists of the electrons 
near the electromagnetic pulse at a; > 0. 



10000 in Fig.|Bfa). Also in Figurel^a), the temporal evo- 
lution of this momentum is shown until tile — 20000 for 
Run A. This temporal increment of the momentum ver- 
ifies that the average electon near the pulse is gradually 
accelerated. In addition, we can see in Fig. Efa) that the 
tail part of the distribution develops into a quasi-power- 
law. At trie — 10000, a second peak of the distribution 
appears at a high value Px, which is consisted of non- 
thermal electrons. In short, some of the surface electrons 
is more accelerated than other surface electrons. This 
multi-peak structure of the momentum distribution may 
be related to the bifurcation of the electron density shown 
in Fig.Etd). 



C. Frequency ratio dependence of the DRPA 

In the previous sections, we summarize the results of 
Run A (with the frequency ratio cjpe/rie — 0.105). Next, 
we examine the frequency ratio dependence of the DRPA. 
The frequency ratios discussed here are cjpe/Oe = 1.0 
(Run B) and Wpe/^^e = 0.01 (Run C). Other fundamental 
simulation parameters are exactly same as those of Run 
A. 

Figure IHl displays the spatial profiles of the same quan- 
tities shown in Fig. [3 in the case of Run B, at different 
times. Figs. IHJa), (b), (c), and (e) mainly show the re- 
gion near the EM pulse at each time, and the values of 
Ez and Jz are zero at tile — 0. The electron density 
rie at tOe = is not shown in Fig. ISfd) since the ini- 
tial density rie/^o = 1 is much larger than the values at 
the other three times. Compared with the results of Run 
A, no apparent EM pulse with large amplitudes can be 
formed at the expansion surface in Figs.|HIa) and |Htb). 
In Fig. ISJ^d), we find that the electrons are distributed 
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FIG. 7: Temporal evolutions of the peak electron momentum 
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FIG. 8: Spatial profiles of the (a) magnetic field By, (b) elec- 
tric field Ey, (c) current density Jz, and (d) electron density 
rie as functions of a; at z = 0, and (e) phase plots of electrons 
in the x — Px space for Run B. The quantities only in the 
a; > region are shown at tQ,e — 0, 500, 1000, and 1500. 



in a broad region, rather than confined in the EM pulse. 
We see no sign of bifurcation of the density pulse like in 
Run A. In Run B, initially the thermal energy of plasma 
particles is much larger than the magnetic field energy 
because of /3e = /?„ = 50. Temporal evolution of each 
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FIG. 9: Temporal evolutions of the peak value of the mag- 
netic field of the pulse for (a) Runs A, B, C, and the case with 
no plasma (pure vacuum EM pulse), and (b) comparison be- 
tween Runs A and D. 



energy component is shown in Fig. ^b). In Run A, the 
field energy is converted into the kinetic energy of sur- 
face particles by the DRPA, thus the total field energy 
decreases as time elapses. In Run B, however, the field 
energy continues to increase as shown in Fig. ^b). This 
increase of the EM field energy is mainly caused by the 
field fluctuation in the regions away from the expansion 
surface. 

We focus on the amplitude of the propagating mag- 
netic field formed at the expansion surface. Temporal 
evolution of the peak value of the magnetic field pulse 
By is shown in Figure El^a). As for Run A, the peak 
value dramatically decreases during < tile < 80 (this 
period corresponds to Phase 1 in Fig.^ and then contin- 
ues to decrease at a constant damping rate (this period 
corresponds to Phases 2 and 3 in Fig.^. We can see the 
similar tendency in the temporal evolution of the peak 
value for Run B during < < 500, but the damping 
rate of the peak value is much larger than that of Run A. 
At trie — 500, the peak value becomes By < 0.3Bq. This 
means that the magnetic field energy in Run B is dis- 
sipated much faster than that in Run A. Although the 
magnetic field energy is rapidly converted into particle 
energy, the total kinetic energy of Run B hardly changes 
since the initial plasma beta is very high. Therefore, 
the evolution can be described as the free-expansion of 
the plasma rather than the expansion accompanied by 
the acceleration of the EM pulse. The electron momen- 
tum peak of Run B, at which the momentum distribu- 
tion takes a maximum, is shown as a function of time in 
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FIG. 10: Spatial profiles of the (a) magnetic field By, (b) elec- 
tric field Ey, (c) current density J^, and (d) electron density 
rie as functions of x at z = 0, and (e) phase plots of electrons 
in the x — Px space for Run C. The quantities only in the 
x>0 region are shown at tQe = 0, 20000, 40000, 60000, and 
80000. 
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FIG. 11: Spatial profiles of the (a) magnetic field By, (b) 
electric field Ey, (c) current density J^, and (d) electron den- 
sity Ue as functions of a; at z = 0, and (e) phase plots of 
electrons in the x — Px space for Run D. The quantities only 
in the a; > region are shown at t^le = 0, 6000, 10000, 14000, 
and 18000. 



Fig. EJb). The momentum remains around — Ibm^c 
even when the amplitude of the magnetic field pulse be- 
comes By ~ 0.3i3o at trie = 500 because the converted 
magnetic field energy is too small to accelerate all the 
plasma particles via the DRPA. Thus, the electrons at 
late times exhibit almost uniform distribution by the free- 
expansion and the density pulse is not formed in Run B, 
as shown in Fig. |HId). We see the evolution of the elec- 
tron momentum distribution for Run B in Fig. EJb) . The 
distribution keeps a Maxwellian form even at late times. 
This means that the plasma is consisted almost only of 
the thermal component. 

Figure EH shows the results of Run C [uupe/^e = 0.01 
and Pe — Pp — 0.005). In this nm the initial magnetic 
field energy is much greater than the thermal energy of 
plasma particles. The converted energy of the EM pulse 
can be effectively used for the particle acceleration [see 
Fig. ^c)]. In Fig. [rUT dl. most of the electrons has been 
accelerated owing to the EM pulse and the electron den- 
sity is much more concentrated at the expanding surface. 
Temporal evolution of the maximum amplitude of the 
pulse magnetic field in Run C is also plotted in Fig.|5fa). 
The period during < tile < 750 represents an initial 
relaxation phase in Run C, which corresponds to Phase 
1 in Fig.^a) of Run A. After that, the peak value stays 
almost constant. By ~ 0.62i?o. For comparison, the re- 
sult of a vacuum EM pulse without plasma is shown in 
Fig. I^a). In this run without plasma, we set up the ini- 
tial magnetic field By to be By{x) = i?o for < 3Aa; 
and By{x) = for \x\ > 3Ax, and the electric field is 
initially zero everywhere. The EM pulse generated by 
the spatial gradient of the initial magnetic field propa- 
gates according to the wave equation in a vacuum. It is 



found that the amplitude of the EM wave is almost un- 
changed after the initial relaxation phase (0 < tfle < 80) 
in Fig. inia). The EM pulse in Run C is very similar to 
the pulse without plasma since the initial plasma beta is 
small, /3 = 0.005, and only a small portion of the total 
magnetic field energy is converted into the plasma kinetic 
energy. 

The evolution of the electron momentum distribution 
for Run C is shown in Fig. I^c). Since most of electrons 
are accelerated by the DRPA, the tail part of the distribu- 
tion exhibits a power-law with sharp slope. In Fig.[2{c), 
the peak momentum of the surface electrons increases 
rapidly as time elapses, compared to Runs A and B with 
higher frequency ratios. The average momentum of elec- 
trons reaches < Pr >~ 1500topC at tflp — 30000. 



D. Dependence on the initial plasma width 

In this section, we examine how the width of the initial 
plasma can affect the DRPA. The widths of the density 
pulse of expanding plasmas tend to scale with the width 
of the initial plasma, as we expect from causality. The 
propagating EM pulses have the same tendency since the 
initial magnetic field has the same width as the plasma 
density. We have carried out Run D, in which the fre- 
quency ratio is ujpe/^le = 0.105 and the initial width of 
plasma is five times that of Run A. 

The temporal evolution of each energy is shown in 
Fig. n^d). The apparent energy conversion from field 
energies into particle energies due to the DRPA can be 
seen until tQ^ — 5000, which is similar to the evolu- 
tion in Fig. nia). After that time, the total EM energy 
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slightly increases and the total particle energies tends to 
decrease. These slight energy changes against the DRPA 
is caused by the local energy changes in the regions away 
from the expansion surface. In those regions, the parti- 
cle acceleration due to the DRPA hardly works efhciently. 
Even though the total particle energies slightly decreases, 
we will show that the plasmas near the propagating EM 
pulse continue to be energized after t^l^ — 5000. 

The electron momentum distributions for Run D are 
displayed in Fig. Eld) at three different times. As time 
elapses, the average momentum of surface electrons in- 
creases, and we recognize a lot of peaks in the high mo- 
mentum part at ifie = 10000. As we will show later in 
Figure ^2 these peaks correspond to a repeated bifur- 
cation of electron density at the expansion surface. In 
Fig-EJa), the results of Run D are also plotted. There 
is no big difference between the peak electron momenta 
of Runs A and D, except in D the acceleration seems 
more delayed due to the prolific early bifurcations. How- 
ever the maximum energy achieved by the most energetic 
particles is higher in D than in A. This is likely caused 
by the larger coherence length for acceleration in D due 
to the wider EM pulse. 

Fig. shows the spatial profiles of the (a) magnetic 
field By and (b) electron density rig, and (c) the phase 
plots X — Px of electrons for Run D at tfle = 6000, 
10000, 14000, and 18000. Only the portion near the ex- 
pansion surface at each time has been displayed in each 
panel. Compared with the electron density of Run A 
in Fig. EJb), the width of the density of the expand- 
ing electrons w becomes larger, e.g., w ~ 2000c/f2e at 
trie = 18000 in Fig. lllf b). We also recognize many more 
bifurcations of the electron density at late times. 

IV. SUMMARY 

We have used two-and-a-half-dimensional large-scale 
particle-in-cell plasma simulations to investigate the 
particle energization of expanding relativistic electron- 
positron plasmas. When magnetized plasmas with high 
temperature expand into a vacuum, the particles at the 
expansion surface are accelerated in the expansion di- 
rection. This energization mechanism is called DRPA 
(diamagnetic relativistic pulse accelerator) Q . 

In this paper, we mainly discuss the long-term evo- 
lution and the initial parameter dependences of DRPA. 
In the case of the electron-positron plasma with the 
frequency ratio Upe/^e — 0.105 and the temperature 
ksTe = ksTp = 5MeV, an electromagnetic (EM) pulse 
with large amplitude is generated and propagates into 
a vacuum. The plasma particles are trapped and ac- 
celerated by the pondermotive force. As time elapses, 
since the magnetic field energy is converted into the ki- 



netic energy of the particles, the amplitude of the EM 
pulse is decreased and the particles near the expansion 
surface are energized. The expanding plasma density is 
concentrated on the vicinity of the expansion surface at 
the early phase of the expansion. And at the late phase, 
we observe the bifurcation of plasma density due to the 
formation of new transverse currents. 

The properties of DRPA strongly depend on the fre- 
quency ratio LUpe/fle- In the case of a high frequency ratio 
{ujpe/ile — 1.0), the EM pulse is rapidly damped and the 
plasma is expanding freely with little acceleration since 
the energy of the plasma is much larger than the mag- 
netic energy. On the other hand, in the case of a low 
frequency ratio {uipe/fle = 0.01), the magnetic energy is 
much higher than the plasma thermal energy. Most of 
the expanding plasma is trapped by the EM pulse and 
the plasma is efRciently energized at late times, achieving 
higher Lorentz factors than in the other cases. 

In addition to the frequency ratio, we have studied the 
effect of varying the initial width of the plasma slab. As 
expected, the width of the asymptotic EM pulse tends to 
be proportional to the initial width of the plasma slab. 
Both the EM pulse and the plasma density are broad- 
ened, and the density pulse develops many more peaks 
at late times due to prolific bifurcations. We also find 
that the maximum energy of the accelerated particles is 
higher in this case due to the longer pulse length. 

Finally, we note the system length in the z direc- 
tion in our two-dimensional simulation domain. As we 
mentioned in Sec. IIII Al the effective Larmor radius of 
relativistic electron-positron plasmas j{c/fle) becomes 
comparable or larger than Lz at late times. Thus, the 
present system length will be insufficient to fully 
cover the physical phenomena scaled by the Larmor 
radius if there exist some instabilities which can de- 
velop along Lz on that scale. However, to our knowl- 
edge, previously observed instabilities in non-relativistic 
simulations 0, such as the lower hybrid drift in- 

stability, do not grow in our electron-positron relativistic 
plasmas. Future advanced simulations will be able to 
give the clear answer to this problem through further 
simulation studies with larger system lengths. 
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